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The magnon exchange mechanism of ferromagnetic super- 
conductivity (FM-superconductivity) was developed to ex- 
plain in a natural way the fact that the superconductivity in 
UGe2, ZrZn-2 and URhGe is confined to the ferromagnetic 
phase. The order parameter is a spin anti-parallel component 
of a spin-1 triplet with zero spin projection. The transverse 
spin fluctuations are pair forming and the longitudinal ones 
are pair breaking. In the present paper, a superconducting 
solution, based on the magnon exchange mechanism, is ob- 
tained which closely matches the experiments with ZrZn^ 
and URhGe. The onset of superconductivity leads to the ap- 
pearance of complicated Fermi surfaces in the spin up and 
spin down momentum distribution functions. Each of them 
consist of two pieces, but they are simple-connected and can 
be made very small by varying the microscopic parameters. 
As a result, it is obtained that the specific heat depends on 
the temperature linearly, at low temperature, and the coeffi- 
cient 7 = ^ is smaller in the superconducting phase than in 
the ferromagnetic one. The absence of a quantum transition 
from ferromagnetism to ferromagnetic superconductivity in a 
weak ferromagnets ZrZn^ and URhGe is explained account- 
ing for the contribution of magnon self-interaction to the spin 
fluctuations' parameters. It is shown that in the presence of 
an external magnetic field the system undergoes a first order 
quantum phase transition. 
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I Introduction 

Very recently ferromagnetic superconductivity (FM- 
superconductivity) has been observed in UGe2 0, 
ZrZn2 Q and URhGe The superconductivity is 

confined to the ferromagnetic phase. Ferromagnetism 
and superconductivity are believed to arise due to the 
same band electrons. The persistence of ferromagnetic 
order within the superconducting phase has been ascer- 
tained by neutron scattering. The specific heat anomaly 
associated with the superconducting transition in these 
materials appears to be absent. 

At ambient pressure UGe2 is an itinerant ferromag- 
net below the Curie temperature Tc — 52K, with low- 
temperature ordered moment of fig — i-^^fJ-s/U. With 
increasing pressure the system passes through two suc- 
cessive quantum phase transition, from ferromagnetism 
to FM-superconductivity at P 10 kbar, and at higher 
pressure Pc '■-^ 16 kbar to paramagnetism [|l],D. At the 
pressure where the superconducting transition tempera- 
ture is a maximum Tgc — 0.8K, the ferromagnetic state 
is still stable with Tc = 32K, and the system undergoes 



a first order metamagnetic transition between two ferro- 
magnetic phases with different ordered moments |^ . The 
specific heat coefficient 7 = ^ increases steeply near 11 
kbar and retains a large and nearly constant value Q . 

The ferromagnets ZrZn2 and URhGe are supercon- 
ducting at ambient pressure with superconducting criti- 
cal temperatures Tgc = 0.29K and Tgc = 0.25K respec- 
tively. ZrZn2 is ferromagnetic below the Curie tem- 
perature Tc = 28. 5K with low-tcmpcrature ordered mo- 
ment of = O.llfiB per formula unit, while for URhGe 
Tc — 9.5K and fig = 0.42/.ts. The low Curie temper- 
atures and small ordered moments indicate that com- 
pounds are close to a ferromagnetic quantum critical 
point. A large jump in the specific heat, at the tem- 
perature where the resistivity becomes zero, is observed 
in URhGe. At low temperature the specific heat coeffi- 
cient 7 is twice smaller than in the ferromagnetic phase, 
materials 

The most popular theory of FM-superconductivity is 
based on the paramagnon exchange mechanism [0,1). 
The order parameters are spin parallel components of the 
spin triplet. The superconductivity in ZrZn2 was pre- 
dicted, but the theory meets many difficulties. In order 
to explain the absence of superconductivity in param- 
agnetic phase it was accounted for the magnon param- 
agnon interaction and proved that the critical temper- 
ature is much higher in the ferromagnetic phase than 
in the paramagnetic one [^j. To the same purpose, the 
Ginzburg-Landau mean-field theory was modified with 
an exchange-type interaction between the magnetic mo- 
ments of triplet-state Cooper pairs and the ferromagnetic 
magnetization density pO| . 

The Fay and Appel (FA) theory predicts that spin up 
and spin down fermions form Cooper pairs, and hence 
the specific heat decreases exponentially at low temper- 
ature. The phenomenological theories |]ll| circumvent 
the problem assuming that only majority spin fermions 
form pairs, and hence the minority spin fermions con- 
tribute the asymptotic of the specific heat. The co- 
efficient 7 = y is twice smaller in the superconduct- 
ing phase, which closely matches the experiments with 
URhGe but does not resemble the experimental re- 
sults for U Ge2 and ZrZn2. The assumption seems to be 
doubtful for systems with very small ordered moment. 

The superconducting critical temperature in (FA) the- 
ory increases when the magnetization decreases and very 
close to the quantum critical point falls down rapidly. It 
has recently been the subject of controversial debate. It is 
obtained in , by means of a more complete Eliashberg 
treatment, that the transition temperature is nonzero at 
the critical point. In fl^ , however, the authors have 
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shown that the reduction of quasiparticle coherence and 
Ufe-time due to spin fluctuations is the pair-breaking pro- 
cess which leads to a rapid reduction of the superconduct- 
ing critical temperature near the quantum critical point. 

Recent studies of polycrystalline samples of UGe2 show 
that T-P phase diagram is very similar to those of single- 
crystal specimens of UGe2 These findings suggest 
that the superconductivity in UGe2 is relatively insensi- 
tive to the presence of impurities and defects which ex- 
cludes the spin parallel pairing. 

Despite of the efforts, the improved theory of param- 
agnon induced superconductivity can not cover the whole 
variety of properties of FM-superconductivity. 

Magnon exchange mechanism of superconductivity has 
been proposed [|l5| to explain in a natural way the fact 
that the superconductivity in UGe2, ZrZn2 and URhGe 
is confined to the ferromagnetic phase. The order param- 
eter is a spin anti-parallel component || -I- || of a spin-1 
triplet (TTiTi + itiii) with zero spin projection. The 
transverse spin fluctuations are pair forming and the lon- 
gitudinal ones are pair breaking. The competition be- 
tween magnons and paramagnons explains the existence 
of the two successive quantum phase transitions in UGe2. 

An itinerant system is considered in which the spin-i 
fermions responsible for the ferromagnetism are the same 
quasiparticles which form the Cooper pairs. Hence, one 
has to consider the equation for the gap as well as the 
equation for the magnetization. Then the system of equa- 
tions for the gap and for the magnetization determines 
the phase where the superconductivity and the ferromag- 
netism coexist. 

In the present paper, a superconducting solution, 
based on the magnon exchange mechanism, is obtained 
which closely matches the experiments with ZrZn2 and 
URhGe. The proposed superconducting solution differs 
from the superconducting solution discussed in in 
two ways. First, the quantum phase transition from fer- 
romagnetism to FM-superconductivity is smooth second 
order phase transition, which resembles the experimental 
results for UGe2- In the present paper, the system passes 
trough a first order quantum phase transition. Second, to 
form a Cooper pair an electron transfers from one Fermi 
surface to the other. As a result, the onset of super- 
conductivity in [ p5[ leads to the appearance of two Fermi 
surfaces in each of the spin up and spin down momentum 
distribution functions. The existence of the two Fermi 
surfaces explains the linear dependence of the specific 
heat at low temperatures as opposed to the exponential 
decrease of the specific heat in the BCS theory. In the fer- 
romagnetic phase the specific heat constant 7 is smaller 
than in the superconducting one, which closely matches 
the experiments with UGe2 § . The onset of supercon- 
ductivity, in the preset paper, leads to the appearance 
of complicated Fermi surfaces in the spin up and spin 
down momentum distribution functions. Each of them 
consist of two pieces, but they are simple-connected and 
can be made very small by varying the microscopic pa- 
rameters. As a result, 7 = y can be made smaller in 



FM-superconducting phase in agreement with URhGe 
experiments H). 

The existence of two Fermi surfaces in each of the spin- 
up and spin-down momentum distribution functions is 
a generic property of a FM-superconductivity with spin 
anti-parallel pairing. An important example is the co- 
existence of ferromagnetism and s-superconductivity in- 
duced by phonons ||lq| . The spin fluctuations are pair 
breaking if the order parameter is spin singlet , and 
superconductivity and ferromagnetism coexist if the spin 
fluctuations are weak, while the magnon-induced super- 
conductivity coexist with ferromagnetism close to quan- 
tum critical point where the spin fluctuations are very 
strong. The present paper and ||l^ describe different 
physical realities which lead to coexistence of ferromag- 
netism and superconductivity. 

The paper is organized as follows. In Sec. II an effec- 
tive Hamiltonian is obtained. In Sec III, the magnon- 
induced superconductivity is discussed. The supercon- 
ducting solution known from |^ is reported to complete 
the investigation. Section IV is devoted to the concluding 
remarks. 



II Effective Hamiltonian 

An itinerant system is considered in which the spin- 
i fermions responsible for the ferromagnetism are the 
same quasiparticles which form the Cooper pairs. The 
effective interaction of quasiparticles Ca-{x){c^ (x)) with 
spin fluctuations has the form 

Hs.f ^ J j (fxc+{x)'^c{x) ■ M{x) (1) 

where the transverse spin fluctuations are described by 
magnons 

Mi{x) + iM2{x) = %/2Ma(f), 

Ml (x) - iM2 {x) = \/2Ma+ [x) (2) 

and the longitudinal spin fluctuations by paramagnons 

Nh{x) - M = ^{x). (3) 

M is zero temperature dimensionless magnetization of 
the system per lattice site. 

The partition function can be written as a path inte- 
gral over the complex functions of the Matsubara time 
T a(T, a:), a"^(r, x), (p{T, x) and Grassmann fimctions 
c^{t,x),c+{t,x) 

ZiP) ^ J D^i{a+,a,ip,4,c,)e-^. (4) 

The action is a sum of free action Sq and part which 
describes the spin-fermion interaction S'int- 



2 











(27r)3 



a^(T, k)a{T, k) + ijj{k)a^{T, k)a{T, k) 

+ ip{T,k)D-^{T,k)ip{T,-k) (5) 

+ c+(t, k)c^{T, k) + ec{k)c^{T, k)ca(T, k) 

where /? is the inverse temperature. The magnon's dis- 
persion is 



where x = (r, x) , Dpm is the paramagnon propagator 
Eq.(^, and Z?,„ is the magnon Green function 



Dra{x) 



uj{k) — pi? 



(6) 



where the spin stiffness constant is proportional to M 
[p — Mpo). The paramagnon propagator in Matsubara 
representation is 



-Dpm(T, fc) = 



duj 



(2^) 



H 



bk^ 



(7) 



The parameter r is the inverse static longitudinal mag- 
netic susceptibility, which measures the deviation from 
quantum critical point. The constants J, po and b are 
phenomenological ones subject to some relations. Fi- 
nally, the spin-up and spin-down fcrmions have the fol- 
lowing dispersion relations: 



/fc2 JM 



Accounting for Eqs(^ and (||) one obtains the following 
expression for the interacting part of the action 



duj d^k e~'''^^+^^^ 
2^(27r)3 iuo + pP 



(12) 



For the purpose of doing analytical calculations it is 
convenient to approximate the four-fermion interaction 
with the static one. To this end I replace the magnon and 
paramagnon propagators Eqs.(|l2|,[7) by static potentials 

2M , , 

(13) 



2MDra{iO,k) -> V^{k) = 

pk^ 

Dpm{uj, k) — > Vpm{k) = 



1 



+ bk^ 

The next step is to represent the spin anti-parallel com- 
posite field c-\ci as a sum of symmetric and antisymmet- 
ric parts. After some algebra one obtains an effective 
four fermion theory which can be written as a sum of 
four terms. Three of them describe the interaction of the 
components of spin-1 composite fields (TTiTi + itjii) 
which have a projection of spin 1,0 and -1 respectively. 
The fourth term describes the interaction of the spin sin- 
glet composite fields ti — it- The Hamiltonians of in- 
teractions are 



d^k, 



C+(fci)c+(fc2) 



8 J Y (2^)' 

C|(fc2 - fc3)c|(fci -I- fca) Vprniks) 



(14) 



S'int = ^ I dr I d-^x 



'2Mct{T, x)ci{t, x)a{T, x) 



+ V2Mc'j^{T,x)ci{T,x)a+{T,x) (9) 

+ [C^iT,x)c]iT,x) - Cj^(T,f)c|(T,f)j ip{T,x) 



The integral Eq.(^) over the bose fields (a,a"'",(^) is 
Gaussian. Integrating them out, using the formula for 
the Gaussian integral [Q, one obtains a representation 
for the partition function in terms of path integral over 
the Grassmann fields. 



Z(/3) = J Dp{c+,c) 



(10) 



The effective fermion action Sef / is a sum of free part 
and resulting four-fermion interaction 5*^4 



S 



J2 



d xid X2 



c^{xi)ci {xi) - c|(a;i)c|(xi) 

(11) 

2 . d'^Xid'^X2C'l^ixi)c^{xi)Dm{xi - X2)c'^ix2)ci{x2). 



Dpm{xi ~ X2) C| {X2)c-\{X2) - C| {X2)ci{x2) 



H 



J2 



ii 



n 



d^k, 



C+{ki)c+{k2) 



(27r)3 

Ci(fc2 - fc3)ci(fci -I- fca) Vpm{kz) 



(15) 



d-'k, 



8 J V (2^)' 



C+(fci)Cj^(fc2) +C+(/ci)c+(/c2) 



(16) 



C]{k2 - k-i)ci{ki -t- /ca) + cx(fc2 - k^)ci{ki + ks) V-ik^) 



d^k,. 



16 J V (27^) 



n 



4(A:i)c+(fc2)-c+(fci)c+(fc2) 



(17) 



CT(fc2 - fc3)c|(fci + k-i) - C|(fc2 - fc3)c|(fci + ks) V+iks) 
where 
V-ik) 



2M 1 

pk^ 



bk^ 



pk^ 



bk^ 
(18) 
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The spin singlet fields' interaction Eq.(|T^ is repulsive 
and does not contribute to the superconductivity 0. 
The spin parallel fields' interactions Eqs.( p^ , p^ are due 
to the exchange of paramagnons and do not contribute 
to the magnon-mediated superconductivity. The relevant 
interaction is that of the ti + it fields Eq.(^6|). It has 
an attracting part due to exchange of magnons and a 
repulsive part due to exchange of paramagnons. 

The effective Hamiltonian of the system is 



^cff — Hq + Hp 



(19) 



where Hq is the Hamiltonian of the free spin up and spin 
down fermions with dispersions Eq.(||). 

Ill Magnon-Induced Superconductivity 

By means of the Hubbard-Stratanovich transformation 
one introduces ti + it composite field and then the 
fermions can be integrated out. The obtained free energy 
is a function of the composite field and the integral over 
the composite field can be performed approximately by 
means of the steepest descend method. To this end one 
sets the first derivative of the free energy with respect to 
composite field equal to zero, this is the gap equation. To 
ensure that the fermions which form Cooper pairs are the 
same as those responsible for spontaneous magnetization, 
one has to consider the equation for the magnetization 
as well. 



M 



1 



< 



> 



(20) 



The system of equations for the gap and for the magneti- 
zation determines the phase where the superconductivity 
and the ferromagnetism coexist. 

The system can be written in terms of Bogoliubov ex- 
citations, which have the following dispersions relations: 



Ei{k) 



JM 



e2(fc) + |A(fc)|^ 



where A(fc) is the gap, and e(fc) — ^ 
temperature the equations take the form 



(21) 

[i. At zero 



M 



1 
2 

J2 



(2^ 

d^k V{p- k) Q{-E2{k)) 



e2(fc) + |A(fc)|2 



A(fc) 



(22) 
(23) 



The gap is an antisymmetric fimction A(— fc) = 
— A(fc), so that the expansion in terms of spherical har- 
monics Yim(S^j:) contains only terms with odd I. I assume 
that the component with I — 1 and m = is nonzero and 
the other ones are zero 



A(fc) Aio(fc)W — cos6'. 

An 



(24) 



Expending the potential V- (k) in terms of Legendre poly- 
nomial Pi one obtains that only the component with I = 1 
contributes the gap equation. The potential Vi (p, fc) has 
the form, 



Viip,k) 



3M 

P 
3M 



P 



'p^ -t- P 



In 



In 



p + fc \ _ l_ 
p — k J pk 

r' + {p + kf 



Ap'^k^ r' + (p- fc)2 



1 

pk 



(25) 



where ^ A, ^ ^ _^ = ^ > i and r' = f « l. 
A straightforward analysis shows that for a fixed p , the 
potential is positive when fc runs an interval around p 
(p — A, p -j- A), where A is approximately independent on 
p. In order to allow for an explicit analytic solution, I 
introduce further simplifying assumptions by neglecting 
the dependence of Aio(fc) on fc (Aio(fc) — Aio(p/) — 
A) and setting Vi(p/,fc) equal to a constant Vi within 
interval {pf — A,nf 4- A) and zero elsewhere. The system 
of equations (E3,E3|) is then reduced to the system 



M = 



dfcfc2 / dt[l - e{-E2{k,t))] 



A 



P/+A 1 

^ f Q{-E2{k,t)) 
327r2 



Pf-A 

where t = cos 9. 



dkk^ I dtt^ 

1 JeHk) + ^t^A^ 



(26) 



(27) 



III.l Solution which satisfies a/-A < JM 



The equation of magnetization (^6|) shows that it 
is convenient to represent the gap in the form A = 
^k(M) JM, where k{M) < 1. Then the equation 



E2{k,t)^0, 

defines the Fermi surfaces. 



(28) 



pf = \^jpj±m^ J^M^ t^A^ , Pf = y/2^m (29) 

The domain between the Fermi surfaces contributes to 
the magnetization M in Eq.(p6|), but it is cut out from the 
domain of integration in the gap equation Eq.(p7[). When 
the magnetization increases, the domain of integration 
in the gap equation decreases. Near the quantum critical 
point the size of the gap is small, and hence the linearized 
gap equation can be considered. Then it is easy to obtain 
the critical value of the magnetization Msc ¥M 
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When the magnetization approaches zero, the domain 
between the Fermi surfaces decreases. One can approxi- 
mate the equation for magnetization Eq.(pq) substituting 
from Eq.(29) in the the difference {pfr ~ (Pj)'^ ^^'^ 
setting = Pf elsewhere. Then, in this approximation, 
the magnetization is hnear in A, namely 



A = J-JkM 



(30) 



where k runs the interval (0, 1), and satisfies the equation 

a^2 



•u\/ 1 — 



mpfj 



(31) 



The Eq.(|3l|) has a solution if mpfJ > 16tt. Substituting 
M from Eq.(^) in Eg. (p7|) , one arrives at an equation 
for the gap. This equation can be solved in a standard 
way and the solution is 



, 16n An fK 

A = \l— i-^exp 

3 m 



oHi^) - 72T7 

2 J'^Vimpf 



(32) 



i 

I{k) = J dtt^ In (f + - nH'' 



Eqs (|^,|l]j32|) are the solution of the system Eqs.(||j27|) 

near the quantum transition to paramagnetism. 

One can write the momentum distribution functions 
{p, t) and (p, t) of the spin- up and spin-down quasi- 

particles in terms of the distribution functions of the Bo- 

goliubov fcrmions 



n^(p, t) 
n^(p, t) 



u^{p,t)ni{p,t) +v^{p, t)n2{p, t) 
u''{p,t){l-ni{p,t))+v''{p, t){l 



(33) 



n2{p,t)) 



where u{p, t) and v(j), t) are the coefficients in the Bogoli- 
ubov transformation. At zero temperature ni(p,t) = 1, 
n2{p,t) = Q{—E2{p,t)), and the Fermi surfaces Eq.(|29|) 
manifest themselves both in the spin-up and spin-down 
momentum distribution functions. The functions are de- 
picted in Fig.l and Fig. 2. 

The two Fermi surfaces explain the mechanism of 
Cooper pairing. In the ferromagnetic phase and 
have different (majority and minority) Fermi surfaces 
(see Fig. 3 and Fig. 4, t = Q graphs). The spin-up elec- 
trons contribute the majority Fermi surface, and spin- 
down electrons contribute the minority Fermi surface. 
When the value of the momentum of the emitted or ab- 
sorbed magnon lies within interval (p/ — A,pf + A) the 
effective potential between spin-up and spin-down elec- 
trons is attracting. Hence, if the Fermi momenta pj: and 

pjp lie within interval {pf — A,pf + A) the interaction be- 
tween spin-up electrons, which contribute the majority 
Fermi surface, and spin-down electrons, which contribute 
the minority Fermi surface, is attracting. As a result, 
spin-up electrons from majority Fermi surface transfer 



to the minority Fermi surface and form spin anti-parallel 
Cooper pairs, while spin-down electrons from minority 
Fermi surface transfer to the majority one and form spin 
anti-parallel Cooper pairs too. As a result, the onset of 
superconductivity is accompanied by the appearance of 
a second Fermi surface in each of the spin-up and spin- 
down momentum distribution fmictions (see Fig. 3 and 
Fig. 4, t = 1 graphs). 




FIG. 1. The zero temperature momentum distribution n, 
for spin up fermions, as a function oi q = ^ and t — cos 6'. 




FIG. 2. The zero temperature momentum distribution n, 
for spin down fermions, as a function of q — ^ and t = cos 6. 
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1 2 ^ 

FIG. 3. The zero temperature momentum distribution n, 
for spin up fermions, as a function of g = — for t — (the 
gap iz zero) and t = ±1 (the gap is maximal). 




1 2 " 

FIG. 4. The zero temperature momentum distribution n, 
for spin down fermions, as a function of q — ^ for t — (the 
gap is zero) and t = ±1 (the gap is maximal). 

The existence of the two Fermi surfaces explains the 
linear dependence of the specific heat at low tempera- 
tures: 



C 
T 



2^2 



{N+{0) + N-{0)) 



(34) 



Here A^^(O) are the density of states on the Fermi sur- 
faces. One ca rewrite the 7 = ^ constant in terms of 
Elliptic Integral of the second kind E{a, x) 



III. 2 Solution which satisfies ,/-^A > JM 



In the present sub-chapter one looks for a solution of 
the system which satisfies 



A > JM 



(36) 



The inequality Eq.(^6|) shows that the gap can not be 
arbitrarily small when the magnetization is finite. Hence 
the system undergoes the quantum phase transition from 
ferromagnetism to FM-superconductivity with a jump. 
Approaching the quantum critical point from the ferro- 
magnetic side, one sets the gap equal to zero in the equa- 
tion for the magnetization (pq) and considers the gap 
equation (|2^) with magnetization as a parameter. It is 
more convenient to consider the free energy as a function 
of the gap for the different values of the parameter M. 
To this purpose I introduce the dimensionless "gap" x 
and the parameters s, A and g 

- — A, s = — JM, A = — , .g = 37 

Then the free energy is a function of x and depends on 
the parameters s, A and g. 

l + A 1 

F{x) = ^ {T{x) - T{0)) = x^+g / dqq^ dt x 

l-A -1 

(s - y/{q^ - 1)2 +t2a;2) {q^ - 1)^ + t^X^ - s)- 



1) 



1) 



(38) 



The dimensionless free energy F{x) is depicted in Fig. 5 
for A — 0.08, g = 20 and three values of the parame- 
ter s, s ~ 0.8, s — 0.69 and Scr — 0.595. As the graph 
shows, for some values of the microscopic parameters A 
and g, and decreasing the parameter s (the magnetiza- 
tion), the system passes trough a first order quantum 
phase transition. The critical values Scr and Xcr satisfy 
? A, 



> 1 in agreement with Eq.(p6| 



mpf 
3k 



(1 + s)-^E{- arcsmK, 

2 s + 1 

(1 — s)^E(- arcsiuK, ) 
^ ' ^2 \s - r 



(35) 



where s — 



JMm 
7^ — 



< 1 and K 



TheEq.(M) 



~ ' ■ ^ y TT JM ' 

shows that in the ferromagnetic phase (A = 0) the spe- 
cific heat constant 7 is smaller than in the supercon- 
ducting one, which closely matches the experiments with 
UGe2. 




0.5 1 1.5 

FIG. 5. The dimensionless free energy F{x) as a function 
of dimensionless gap a:. A = 0.08, g — 20, si = 0.8(upper 
line), S2 = 0.69(middle line) and s^r ~ 0.595(lower line). 
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Varying the microscopic parameters beyond the critical 
values, one has to solve the system of equations ( ^^ . 
One represents again the gap in the form 



(39) 



A = \ j-K{M)JM 



but now k{M) > 1. Then the equation E2{k,t) — 0, 
which defines the Fermi surface, has no solution if — 1 < 



t < 



when 



k(M) 



and 



k(M) 



< t < 1, and has two solutions 



3 

— ! 

TT 



(40) 



k{M) 

The solutions 



< t < 



k{M) ■ 



determine the two pieces of the 



k(M) ' 



SO 



Fermi surface. They stick together at t — ± 
that the Fermi surface is simple connected. The domain 
between pieces contributes to the magnetization Af in 
Eq. ( p6| ) , but it is cut out from the domain of integration 
in the gap equation Eq. (p7|) . The Fermi surface manifests 
itself both in the spin-up and spin-down momentum dis- 
tribution functions. The functions are depicted in Fig. 6 
and Fig. 7. 

When the magnetization approaches zero, one can ap- 
proximate the equation for magnetization Eq. (Eq) substi- 
tuting from Eq.(^ in the the difference (pjr — (pj)^ 

and setting = Pf elsewhere. Then, in this approxima- 
tion, the magnetization is linear in A, namely 



A = J-JkM 



(41) 



where n = H^mJ. jg the small magnetization limit of 
k{M) The Eq.(^) is a solution if mpfj > IGtt (see 
Eq.(|6|)). Substituting M from Eq.(|l|) in Eq.(|7|), one 
arrives at an equation for the gap. This equation can be 
solved in a standard way and the solution is 



16tt pfA 



exp 



(42) 



Eqs (^^2|) are the solution of the system near the quan- 
tum transition to paramagnetism. The second derivative 
of the free energy Eq. ( |38| ) with respect to the gap is pos- 
itive when "t^/'^ > (^rs^)^, hence the state where the 
superconductivity and the ferromagnetism coexist is sta- 
ble. 

The existence of the Fermi surface explains the linear 
dependence of the specific heat at low temperature: 



C 
T 



NiO) 



(43) 



Here A''(0) is the density of states on the Fermi surface. 
One can rewrite the 7 = y constant in terms of Elliptic 
Integral of the second kind E{a, x) 



mpf 
3k(M) 



(l + ^)^i^(4,^)+ 



> 1 „,7r 2s . 



(44) 



where s < 1 (see Eq.(| 




FIG. 6. The zero temperature momentum distribution n, 
for spin up fermions, as a function of g = ^ and t = cos 6. 




FIG. 7. The zero temperature momentum distribution n, 
for spin down fermions, as a function of q — ^ and t — cos 6. 

Eq.(Q) shows that for k{M) just above one the spe- 
cific heat constant 7 is smaller in ferromagnetic phase, 
while for k{M) » 1 it is smaller in FM-superconducting 
phase. The result closely matches the experiments with 
ZrZn2 and U RhGe respectively. 

The solutions Eqs . ( pO |4l| ) show that magnetization and 
superconductivity disappear simultaneously. It results 
from the equation of magnetization, which in turn is 
added to ensure that the fermions which form Cooper 
pairs are the same as those responsible for spontaneous 
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magnetization. Hence, the fundamental assumption that 
superconductivity and ferromagnetism are caused by the 
same electrons leads to the experimentally observable 
fact that the quantum phase transition is a transition 
to paramagnetic phase without superconductivity. 

An important experimental fact is that ZrZni and 
URhGe are superconductors at ambient pressure as op- 
posed to the existence of a quantum phase transition 
in UGe2. To comprehend this difference one considers 
the potential (|25|). The quantum phase transition re- 
sults from the existence of a momentum cutoff A, above 
which the potential is repulsive. In turn, the cutoff ex- 
cistence follows from the relation (3 — > 1, which 
is true when the spin-wave approximation expression for 
the spin stiffness constant p — Mpo is used. The spin 
wave approximation correctly describes systems with a 
large magnetization, for example UGe2- But in order to 
study systems with small magnetization, one has to ac- 
count for the magnon-magnon interaction which changes 
the small magnetization asymptotic of p, p = 7\f^+"poj 
where a > 0. Then for a small M (] < 1, and the 
potential is attractive for all momenta. Hence for sys- 
tems which, at ambient pressure, are close to quantum 
critical point, as ZrZn^ and URhGe, the magnon self- 
interaction renormalizes the spin fluctuations parameters 
so that the magnons dominate the pair formation and 
quantum phase transition can not be observed. But if 
one applies an external magnetic field, the magnon opens 
a gap proportional to the magnetic field. Increasing the 
magnetic field the paramagnon domination leads to first 
order quantum phase transition. 



IV Conclusions 

The proposed model of ferromagnetic superconductiv- 
ity differs from the models discussed in [fzt^ljgljiotl in many 
aspects. First, the superconductivity is due to the ex- 
change of magnons, and the model describes in an unified 
way the superconductivity in UGe2, ZrZn^ and URhGe. 
Second, the paramagnons have pair-breaking effect. So, 
the understanding the mechanism of paramagnon sup- 
pression is crucial in the search for the ferromagnetic 
superconductivity with higher critical temperature. For 
example, one can build such a bilayer compound that the 
spins in the two layers are oriented in two non-coUinear 
directions, and the net ferromagnetic moment is nonzero. 
The paramagnon in this phase is totally suppressed and 
the low lying excitations consist of magnons and addi- 
tional spin wave modes with linear dispersion e(fc) ^ k 
pof . If the new spin- waves are pair breaking, their effect 
is weaker than those of the paramagnons, and hence the 
superconducting critical temperature should be higher. 
Third, the order parameter is a spin antiparallel com- 
ponent of a spin triplet with zero spin projection. The 
existence of two Fermi surfaces in each of the spin-up 
and spin-down momentum distribution functions leads 
to a linear temperature dependence of the specific heat 



at low temperature. 

The proposed model of magnon-induced superconduc- 
tivity does not contain the relativistic effects, namely 
spin-orbital coupling which is present in UGe2- The 
resulting magneto-crystalline anisotropy will modify the 
spin-wave excitation and will add a gap in the magnon 
spectrum, which changes the potential Eq. (psf) . The 
physical consequence of the change is that the supercon- 
ductivity disappears before the quantum phase transi- 
tion from ferromagnetism to paramagnetism (see ) . 
The distance between these two points depends on the 
anisotropy. 
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